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I. Introduction 

In Ref. 1, it was presented a method, primarily conceived as an emergency backup system, that addressed 
the problem of a space capsule that needed to execute a safe atmospheric entry from an arbitrary initial 
attitude and angular rate in the absence of nominal control capability. The proposed concept permits the 
arrest of a tumbling motion, orientation to the heat shield forward position and the attainment of a ballistic 
roll rate of a rigid spacecraft with the use of control in one axis only. To show the feasibility of such concept, 
the technique of single input single output (SISO) feedback linearization using the Lie derivative method 1 2 
was employed and the problem was solved for different number of jets and for different configurations of the 
inertia matrix: the axisymmetric inertia matrix {I xx > I yy = Izz ), a partially complete inertia matrix with 
I xx > Iyy > Izz,Ixz 7 ^ 0 and a realistic complete inertia matrix with I xx > I yy > I zz ,Iij yf 0. The closed- 
loop stability of the proposed non-linear control on the total angle of attack, 0, was analyzed through the zero 
dynamics of the internal dynamics 3 for the case where the inertia matrix is axisymmetric ( I xx > I yy = I zz ). 
This note focuses on the problem of the diagonal non-axisymmetric inertia matrix ( I xx > I yy > I zz ), which 
is half way between the axisymmetric and the partially complete inertia matrices. In this note, the control 
law for this type of inertia matrix will be determined and its closed-loop stability will be analyzed using the 
same methods that were used in Ref. 1. In particular, it will be proven that the control system is stable in 
closed-loop when the actuators only provide a roll torque. 

II. System Equations 

The equations describing the control problem are those of a rotating rigid body with extra terms describing 
the effect of the control torques. They therefore consist of kinematic equations relating the angular position 
with the angular velocity, and dynamic equations describing the evolution of angular velocity. 

A. Dynamic Equations 

The dynamics of the rotational motion of a rigid body are described by the Euler’s equations. Let / be 
the inertia matrix of the spacecraft, and let to denote the angular velocity vector with components along a 
body-fixed reference frame located at the center of gravity (c.g.) and aligned along the principal axes of the 
spacecraft. The dynamic equations are presented in Eq. (1), where ku denotes the control torque vector. 

( 0 to z —u), 

-u z 0 u> x 

u) y a) x 0 

B. Kinematic Equations 

The kinematic equations relate the components of the angular velocity vector with the rates of a set of 
parameters that describe the relative orientation of two reference frames: an inertial frame and the body- 
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fixed frame. In the formulation of this problem, the inertial frame OXy ,Yy ,Zy is associated to the velocity 
vector. It is defined with its origin located at the c.g. of the spacecraft, with the Xy axis pointed along the 
velocity vector, the Yy axis pointed along the angular momentum of the trajectory, and Zy axis pointed 
along Yy x Xy. The body-fixed frame Ox,y,z is defined by the body axes. In this case, when the rotation 
angles are zero, the x axis goes along the longitudinal axis, aligned with Xy, being positive toward the heat 
shield, y is aligned with —Yy and z with Zy. For convenience, we will consider a 1-3-1 sequence of rotations 
which is a type 2 Euler sequence. The rotations involved are ip about Xy, 9 about Z' v and ip about x, where 
Zy is the resulting Zy after the first rotation. This selection results in the kinematic equations, Eq. (2) 

p = (— ojy cos ip + u) z sin ip) esc 9 

9 = w y sin ip + co z cos ip (2) 

ip = lo x + ( u) y cos ip — oj z sin ip) cot 9 


C. Equations of Motion 

Combining Eq. (1) and Eq. (2) yield the attitude control problem. We can see that 9, the variable that we 
want to control first, only depends on u y , ui z and ip in the kinematic equations, therefore the equation for 
ip can be ignored. This makes sense since 9, as well as its angular components, angle of attack and sideslip 
angle, are invariant to a rotation around the velocity vector. Hence the equations of motion can be reduced 
to 


Co = / - 1 w x /d; + I~ 1 hu 

9 = co y sin ip + co z cos ip (3) 

ip = lo x + ( u)y cos ip — oo z sin ip) cot 9 

III. Nonlinear SISO System 

The system in Eq. (3) can be written in the form of a SISO system 

X = ( L0 X LOy LO z 9 Ip ) T 

( /" 1 w x /w \ 

f(x) = coy sin ip + co z cos ip (4) 

y lo x + (u> y cos ip — lo z sin ip) cot 9 J 
g(x) = ( I~ x k 0 0 ) T , h(x) = 9 - 9 d 

with x £ SR”, where n is the system order (5 in our case), with / and g smooth vector fields on 3?" and h 
a smooth nonlinear function, and where 9d is the desired, or targeted 9. A large class of SISO nonlinear 
systems can be made to have linear input-output behavior through a choice of nonlinear state feedback 
control law that is given with generality by 

where Lfh(x) and L g h( x) stand for the Lie derivatives of h with respect to / and g respectively, 7 is the 
relative degree of the system and the input v is a scalar that results from the product of the feedback gain 
vector times the error vector. This control law yields the 7-th order linear system from input v to input 
y : '-jpp = v, where 7 is the smallest integer for which LgL^h (x) = 0, i = 0, ...,7 — 2. It is easier to understand 
7 as the smallest integer for which the control signal appears in cTy/dt 7 for the first time. 

By inspection of Eq. (4) we can see that the relative degree 7 will depend on the specific form of the 
product I~ 1 k. Therefore, both, the inertia matrix and the geometric layout of the jets are the two factors 
that affect the relative degree of the system. Note that 7 can only take the values 2 or 3: it cannot be 1 
since the first derivative of the output, 9 , does never contain the control; it is 2 if the vector I~ 1 k has either 
its 2 nd or 3 rd component different than zero; and 7 would be 3 if the vector I~ 1 k only had its 1 st component 
different than zero. A 7 higher than 3 is not possible. 


x = f(x) +g{x)u 
y = h(x) 
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The relative degree, and hence the product has important implications in the stability character- 

istics of the system under control since the dimension of its internal dynamics is given precisely by n — 7. 
When 7 is defined and smaller than n, the nonlinear system in Eq. ( 4 ) can be transformed into a so called 
normal form, which shall allow us to take a formal look at the stability of the system through the notions of 
internal dynamics and zero dynamics. The normal states become 4>(;r) = ( ^ • • • £ 7 rji • • • rj n - 7 ) T 

and the normal form of the system can be written as 


Ci = L f h(x) = C 2 
6 = L}h(x) = C 3 


Vi = Lfi 71 = qi{£,r}) 
1)2 = LfV 2 = 52(C) »?) 


C 7 = L}h(x) + LgLj-'h^u ~ Qn ~ l{ ^ V) 

where the vector field 77 is a solution of the set of partial differential equations 


( 5 ) 


^Vj(x)g{x) = 0 1 < j < n — 7 (6) 

To show that the nonlinear system in Eq. ( 4 ) can indeed be transformed into the normal form of Eq. 
( 5 ), we have to show that we can construct a local diffeomorphism 4 >(a;) such that Eq. ( 5 ) is verified. To 
show that < l>(a;) is a diffeomorphism, it suffices to show that its Jacobian d$/dx is invertible, i.e., that the 
gradients V£i and V?7j are linearly independent (span(d&( x)/dx) € 5ft 5 ). 

The internal dynamics associated with the input-output linearization correspond to the last n — 7 equa- 
tions of the normal form and they constitute the unobservable part. However, the control must account for 
the stability of the whole dynamics and, therefore, the system will be stable if the internal dynamics remain 
bounded. To assess the stability of the internal dynamics we will use the zero dynamics of the system, which 
is defined as the internal dynamics of the system when its output is kept at zero by the input. As mentioned 
in Ref. 3 , the zero dynamics is an intrinsic feature of a nonlinear system that does not depend on the choice 
of the control law or the desired trajectories. 

The constraint that the output h is identically zero in zero dynamics implies that all of its time derivatives 
are zero. Thus, the corresponding internal dynamics of the system, or zero dynamics, describes a motion 
restricted to the n — 7 dimensional manifold defined by £ = 0 . Also, in order for the system to be in zero 
dynamics, the input control u must be such that h stays at zero. This means that in zero dynamics the 
normal form can be written as (£1 = 0 , . . . , £ 7 = 0, 1)1 = <7i(0, 77), . . . , ?)„_ 7 = q„_ 7 (0, 77), h = 0) 

The analysis of the zero dynamics of the internal dynamics will be carried out for the control on 0 in 
the axis-symmetric inertia case only. The resulting expressions in the intermediate and realistic inertias are 
practically intractable. Nonetheless, we will see that useful conclusions can be from the axis-symmetric case 
can be extended to the other more realistic cases. 


IV. Closed-loop Stability with Roll Torque 

In the case in which we have a diagonal non-axisymmetric matrix of inertia and in which only a roll 
torque is provided, the equations of motion become 


( 6 j x \ 


/ au> x u) z \ 


/ I - 1 

XX 

LUy 


buj x u> z 


0 

U) z 

= 

ClVx&y 

+ 

0 
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Lo y sin 1 p + uj z cos ip 


0 

\ ip ) 


\ lo x + (uiy cos ip — co z sin 7 p) cot 9 / 


V 0 


( 7 ) 


where CL — (tyy Izz) / 1 xx , b — ( 1 -zz Ixx) / lyy and C — (Cx tyy) / '^zz- 

The relative degree of the system, 7, is 3 and, thus, the control law will be given by u = (v—IxjO) / ( L g L 2 j. d ), 
which is shown in Appendix A. 

The external part of the normal system is given by 
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£i = $1 = h(x) = 0-0 d 
£2 = $2 = Lfh(x) = 6 
£3 = $3 = Ljh(x) = 6 

where Lfh(x) stands for the Lie derivative of h with respect to /. 

The internal part of the normal system will have 2 components and must satisfy \7r/j(x)g(x) = 0 

^2 ! „( T ) — dm t- 1 — n 

d x y\ x ) — dui^ xx ~ u 

dm n ( T ) - drnj-l - n 

d x y\ x ) — dujx 1 ^ ~ u 

The solution of the internal system must result in = ( £1 £ 2 £3 Vi V 2 ) T being a diffeomor- 

phism. One solution that satisfies this requirement is given by 771 = to y , 772 = to,. Therefore, the internal 
system becomes 


m =u y = buj x u z , g ^ 

772 = X> Z = CU) X U)y 

In the axisymmetric case, it could be demonstrated that u> x was constant in zero dynamics. However, 
that is not necessarily the case now. To check whether the system in Eq. 8 is stable, u> y is divided by oj z , 
resulting in 


LOy_ 

LU Z 


Operating on Eq. 9 and integrating results in 


bui z 

CbJy 


(9) 


(jJydiOy / LU Z diU Z = 0 


w; - -wi = c 

y c 


(10) 


where C is the integration constant which includes the factor 1/2 coming from the integration. Eq. 10 can 
be rewritten as 


LO 


V 

C 


to. 


C{—c/b) 


= 1 


(11) 


which is the standard equation of an ellipse. Note that in Eq. 11 the fraction c/b is negative because 

x xx ^ y u ^ zz • 

This mathematical development proofs that the internal system in Eq. 8 is bounded and thus, that the 
closed-loop control system can be said to be stable. 


A. Test Case 

In the case of an axisymmetric matrix of inertia the transversal rate lot = tor 0 = ( w yo + w zo) 1 ^ 2 was constant 
throughout the control process. Now, this condition does not hold. Nevertheless, the settling time will be 
determined following the logic that was used in the axisymmetric case presented in Ref. 1. A settling time 
t s = (2n/u>To)5 is selected, where S is a factor that represents some percentage of added or subtracted time. 
In this case, <5 = 1.25, which represents a 25% time addition. Nevertheless, for all practical purposes, the 
important time is the one that marks the moment where 6 becomes confined below the side wall angle, which 
is always smaller than t s . 

For the design of the linear controller, the feedback gain vector K is calculated using a pole placement 
technique. It is selected such that a desired set of closed-loop pole locations p is matched. The location of 
the poles was chosen such that there was no overshoot in 9 , at least when the control signal is not saturated. 
The location of the poles will vary with lot 0 in the following manner: 

P ={- 1 -5 -7)5 /t s (12) 

Figure 1 presents the results for the case in which the initial rates are at the edge of a typical rate 
deadband limits. The moments of inertia are given by I xx = 5900 kg m 2 , I yy = 5100 kg m 2 and I~ z = 4700 
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kg m 2 . This case has the following initial conditions: u x o = 2 deg/sec, u y o = -2 deg/sec, w z q = -2 deg/sec, 
do = 150 deg, tjjo = 120 deg, 0d = 9 deg, t B = 159.1 sec. Figure 1 is divided in two parts. The left side 
presents the locus of the positive direction of the spacecraft’s longitudinal axis (+x axis) in the projection 
of the two dimensional attitude sphere in terms of the angles of attack and sideslip. Each of the contours 
represents an iso-0 trajectory spaced every ten degrees around the velocity vector, located at the center of 
the diagram. A spacecraft exactly flying heat shield forward would show its locus right at (0,0), whereas a 
spacecraft exactly flying apex forward would show its locus at (0,±180). The ring in the center represents 
a typical side wall angle of 32.5 deg, therefore, as a first approximation, to have the capsule protected from 
the heat of the entry means that 6 should remain confined within that ring. The right side of the chart 
shows the time histories of the control signal, the total angle of attack and the roll rate. The control signal 
is expressed as a percentage of the maximum acceleration that the system can generate; the plot showing 6 
also indicates the side wall angle as a dashed horizontal line and the settling time t s : and the roll rate plot 
also shows the closest ballistic roll rate value | pt, | or — | pt, |. The case considered in this note is that of 
a ballistic capsule because there are no products of inertia and, thus, the center of gravity is contained in 
the axis of symmetry of the spacecraft. In any case, the roll rate is shown because it is a byproduct of this 
type of control acting on the roll axis. Should the resulting roll rate be so high that it violates the medical 
constraints on angular rates or a possible constraint on the deploy of the landing system, the roll rate could 
be controlled in the presence of atmosphere as it is explained in chapter VII in Ref 1. 




Figure 1. Test case results. 


It can be noted in Figure 1 that the control signal becomes periodic in steady state. Nevertheless, the 
control can be turned off after t = t s , a point where the error signal is close to zero. That condition ensures 
the total angle of attack confinement and a ballistic entry. Figure 2 shows the results when the control signal 
is turned off at t = 175.0 sec. 


V. Conclusion 

It has been shown that for the diagonal non-axisymmetric matrix of inertia with relative degree 3, there 
exists a feedback control law that can transfer the system to a desired total angle of attack with 0 = 0 = 0 
while the rest of the state variables remain bounded. A feedback control law based on Feedback Linearization 
that satisfies the above condition has been formulated. The controller’s gains were designed to depend on 
the initial conditions and they were generated following a standard pole placement methodology. 

Also, for the diagonal non-axisymmetric inertia matrix, the closed loop stability of the system when 
controlling 0 with roll torque was demonstrated to be globally stable. 
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Max Control Accel. = 5m/s 2 






Figure 2. Test case results. Control turned off at t = 175.0 sec. 


Appendix A - Derivation of the Control Law 


The control law u is given by 


where 


LfO = (dO/dx)f(x) 

L“jd = ( dLf6/dx)f(x ) 

L 3 fO = (dLrjO /dx)f{x) 


L g L)6 = ( dL 2 j:d/dx)g(x ) 


u L 3 fQ)/(L g L 2 fO ) 


(13) 


u>y sin ip + to z cos ip 

ui x [wy cos ip(c + 1) + u z sin tp(b — 1)] + (u> y cos ip + w z sin ■ ip ) 2 cot 9 
[u) y cos ip(c + 1) + u> z sin ip(b — 1 )\Cj x 
[lo x cos ip(c + 1) + 2 (ip — uj x ) cos ip]Cj y 
[lo x sin ip{b — 1) — 2(ip — u> x ) sin ip\io z 


<p 2 9 

[lu z cos ip(b — 1) — uj y sin ip(c + 1)] — 2 (ip 
I xx l^y cos + 1) + u z sin ip(b— 1)] 
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